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Analytical Study of Twin-Jet Shielding 
I SUM ARY 

This report summarizes the development of the analytical model of twin- jet 
shielding. The model consists of a point noise source impinging on a cylinder 
of heated flow in which the temperature and flow velocity are uniform across the 
cross-section. 

In the formulation of the model, the wave equations are written for the re- 
gions outside the flow and within the flow cylinder. The solutions to the wave 
equations are matched at the jet boundary under the conditions of continuity of 
pressure and continuity of the vortex sheet. The solution reduces to an indef- 
inite integral involving Bessel functions. The integral is solved approximately 
using a saddle point method. 

The resulting mod<d is analysed in order to identify the mechanisms of trans- 
mission and diffraction which control shielding in the shadow of the shielding jet. 
The impact that variations in jet operating parameters has on shielding i;? inves- 
tigated. It is found that in the zone of the shadow region in which transmission 
dominates, shielding is relatively insensitive to variations of such parameters 
as Mach Number and spacing ratio. Comparison of the three-dimensional model to a 
previously obtained two-dimensional model shows good agreement. 

In the zone in which diffraction dominates, shielding is more sensitive to 
variations in Mach Number, jet temperature and spacing ratio. 

The model is found to compare to experimental shielding results > The major 
discrepancy between the model and experiment is felt to arise from the difference 
between the single point noise source of the model and the distributed source of 
the jet. Refinement of this noise source to a more realistic representation of 
the real jet is one of the goals of the continued effort in this project. 


I INTRODUCTION 


One of the drawbacks of the growing dependence on air travel is the increased 
impact of aircraft noise. Assessment and reduction of this impact requires iden- 
tification of aircraft generated noise levels. To this end, the Noise Technology 
BrancJi of MA&VAircraft Noise Reduction Division is developing and refining an 
aircraft noise prediction computer program. Noise estimation includes consideration 
not only of noise sources on the aircraft, but also of the propagation path between 
source and receiver, One of the numeious factors affecting the noise transmission 
path is shielding of one jet by another. The shielding jet, because of the high 
temperature and flow speed with respect to the immediate surroundings, acts as a 
partial barrier between the source and the receiver. 'The resultant noise reduction 
not only affects the overall aircraft noise level, but also indicates the possi- 
bility of jet engine installation as a means of aircraft noise control. 

It is the pui*pose of this project to develop an analytical method to estimate 
the shielding of one jet by an adjacent jet in a twin jet configuration. 

The problem of reflection and transmission of sound by a moving medium has 
been addressed assuming a plane wave incident on a plane interface ( 1,2, 3, 4 )• 

Ray tracing techniques have been applied to two-dimensional jets (5) and cold 
jets C <j). 

In this study, the three-dimensional case is considered. The model developed 
consists of the sound field emmited from a stationary, discrete frequency point 
source, which impinges on a cylinder of locally parallel flow. The temperature 
and velocity profiles are uniform in the jet at any location down stream of the 
nozzle . While this model is idealized, it is felt to incorporate the basic 
elements essential for a realistic representation; not only of the source, but also 
Of the shielding jet. 
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The wave equations are written for the acoustic velocity potential in the 
region outside the jet and in the region within the jet, The equations incorpo- 
rate the convective effect due to local mean flow, The velocity potentials in 
the two regions are matched at the flow boundary by the conditions of pressure 
continuity and continuity of the vortex sheet. 

Hie solution for the acoustic pressure in the far field, which is derived 
from the velocity potential, is normalized by the incident acoustic pressure. 
This incident sound pressure is the sound pressure at the same location in the 
far field in the absence of the jet. This normalized sound pressure is related 
to the directivity, 

The model is used to evaluate the effect on the far field noise radiation 
due to variations of such operating parameters as shielding jet flow speed and 
temperature, and the spacing between the source and the shielding jet. 


Ill MODEL DEVELOPMENT 
1, Formulation of the Model 

The derivation of the model was shown in the first progress report ( 7 ). 

A swinery of that derivation is shown below. 

The mechanisms by which shielding occurs are reflection and refraction of 
sound at the boundary between the jet and the surrounding air and by diffraction 
around the jet. 

The noise source is modelled by a stationary, discrete frequency point 
source located at (r Q , q q , 6) . The shielding jet is a cylinder of radius a, and 

is infinite in extent along the z-axis . The temperature and flow velocity are 

* 

uniform across the cylinder cross -section. The model is illustrated in Figure 1. 
The expression for acoustic velocity potential is written for two regions; region 
I is outside the jet, region II is within the jet. 

In region l (outside the flow-incident upon the shielding jet) 

V 2 * - • Q 0 *’ lwt $(r-r Q ) 6(e-9 0 ) 6(z) * la) 


In region I (outside the flow-reflected from the shielding Jet) 


- T * 


c. 


tt 


■ 0 


In region II (inside the flow) 


•v 2 ^ - m2 * 2Z 


2M * m ±. 

Cl 9 Zt c 2 9 tt 
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Where: 

(r Q , 9 Q , o,) - location of point source 
(r, 9, z) - location of the receiver 
u> - source frequency 
Q o - source strength 


lb) 


2 ) 


* 


c - sound speed 

M - mach nunber - (jet flow speed/ci) 

* 

Note: The subscript o refers to con4itions outside the flow (ambient) , and 1 
refers to conditions within the heated jet. 

The boundary conditions at the interface between the ambient air and the 
jet are: 

1) Pressure continuity 

(P) 0 ■ (P ) x at r ■ a 
or- 

" ' p i^t + at r ■ a 3) 


2) Continuity of the vortex sheet (1,2 ).This condition states that the 
displacement of the medium is continuous and symmetrical at the boundary; 
r ■ a. Denoting this displacement by n ■ (2 ,t) , then: 


or- 




at r ■ a 


(n t ) 0 ■ (n t + v n s ) i, at r ■ a 


4) 


Time is eliminated from equations 1 and 2 ‘by assuming; 

♦(r, 8, 2 , t) » <J»(r, e, z)e’ lut 


The problem is reduced to a two-dimensional formulation by the Fourier 

*• 

ipe 


transform: 
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z dz 
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with corresponding inverse: 



Solution of the transformed equations , inclusion of the boundary conditions 
and inverse transformations yields the equation for the acoustic velocity 
potential in the far field (r»a) : 


•iQo 1 
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2. Solution, of the ftodel 

An approximate solution of the integral in equation 5 is obtained using 
the Method of Stationary Phase (g). By this method , the solution of the integral 
I z , of the form: 


is: 


I, ■ 


l z ■ 


| hr 


l * |h M MJ 


as z + • 


g (« 0 ) 
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where: 

1\ a Q solves h' (a) ■ o 

2) the sign in the exponential term goes as the sign of h" (t 0 ) 


In OTder to solve equation .5 in the manner prescribed by equation 6 , 
the following transformations are made. 

Let: 


1 / 


k z 2 
1 -C— ) 


then: 


1) dk, 
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the integral in equation 5 : 


V* ' 

■f- v 


,W V k o, k ‘) ■** 


becomes : 






ik Q z fl-a" 


W> da 8) 


In the acoustic far field, r»l. Assuming that k 2 «k Q , the argument of the 
Hankel function becomes large: 


k Q ra» 1 

and the Hankel function is approximated by: 




s i/^ 


i(k rot - V ( m + W 


V a 


9) 


from figure 1, it is seen that: 

r ■ R cos 9 and z ■ R sin 9 

where: 

R ■ distance from origin of the coordinate system to the receiver 
9 • angle between vector R and the projection of the vector R on the 
x - y plane. 

The integral is then: 



comparing equation 10 to equation 6 : 


i *k 0 R 


h(«) ■ Vl-a^sin S * a wJ 8 

•» ■ 


a is evaluated from: 
o 


h* (aj * cos 8 - ==- 
0 


sin 8 - o 


thus: 


a ■ cos 8 
0 


then: 

and: 


h(* t j - 1 


h" (« 0 ) 


sin 2 8 


and: 



cos 8 ' 
sin 2 8 


F m (cos 8 ) 


The integral, as k Q R » 1, becomes: 
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From which the expression for the far field acoustic velocity potential 


is evaluated: 

Oe"^ t • -imir ik R 
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Where: 


F a (co* () ■ J m (k Q r 0 cos B) 
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where: 


T i ■ ^ - M sin s 


T 2 - [ ( — - M sin 8) 2 - sin 2 a] 


The total acoustic velocity potential in the far field consists of two 
1. The velocity potential of the incident wave: 


0 e’i^t, 
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2. The velocity potential of the scattered wave: 

-iwt 
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terms in equation 14) within brackets 
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16 ) 


As a check; when r Q ■ o (source is located at the origin ) , the velocity 
potential of the incident wave reduces to: 


sine*; 


1 m ■ o 


This is the expression for the velocity potential in the far field for 
radiation from a point source of strength, Q 0 , located at the origin of the co- 
ordinate system. 
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IV RESULTS 


The total sound pressure is evaluated from the acoustic velocity potential 
by: 

P T " " p o, ^ $T 

The incident sound pressure is evaluated from equation 15: 

3*. idiQ oo 1^21 ik R 

P in ' - p o T~ ■ * -^r- p o r e m cos ra < 0 -®) e ' e ° Wo cos 6 > 17 > 

444 at 4nR m*o 

The normalized sound pressure is formed from the magnitude of the ratio of 
the total sound pressure to the incident sound pressure. This normalized sound 
pressure is a measure of the influence of the shielding jet on the noise source 
directionality. A value of normalized sound pressure less than one indicates noise 
reduction; and a value of normalized sound pressure greater than one indicates 
amplification. 

The normalized sound pressure is evaluated in the shielded zone of the jet. 

Thus, 0-0 Q = it and o < 8 < tt/z . For the purposes of comparison to other data, 
the coordinates are shifted to a set centered in the source, and illustrated in 
figure 2. I« the nomenclature adopted in the figure, iji = o° is on the z-axis of 
the noise source, and i|/ n * 90° is on the side of the jet directly opposite the 
source. 

1. Identification of Shielding Mechanisms 

The mechanisms which control the noise in the shadow of the shielding jet in 
the downstream region are transmission through the jet and diffraction around the 
jet. 
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In order to differentiate among the mechanisms of shielding, the normalized 
sound pressure is evaluated at low frequencies under the following operating 
conditions, 

a) Vvj * Jet velocity - 1552 feet/sec 

b) Tj ■ Jet tmeperature ■ 1238°R 

c) T q ■ Ambient temperature ■ 530°R 

d) S/D ■ Spacing parameter - ratio of distance from source to center of jet 

and jet diameter. * 2,667 


The contours of normalized pressure for 10° < < 90° are plotted against 

non-dimensional frequency parameter in figure 3. The non-dimensional frequency 

k.a is varied from 0.1 to 1.0 

o 

For angles, # , within the range: 

90° > i j» n > 70° 

transmission of noise through the jet is dominant. The normalized sound pressure 
is close to 1.0 for a relatively large range of frequencies, indicating that the 
noise transmits easily through the jet. The normalized sound pressure curves shift 
toward lower frequency as the angle is decreased, which indicates that as the angle 
of incidence becomes more oblique, more of the incident sound is entrained within 
the jet and and convected downstream, 

The cut-off angle: 


♦net * C0S 


-1 


c./c- 

(A~±) 

1+M 


where : 


o q ■ sound speed at the ambient temperature, T Q 
c. » sound sieed at jet temperature, T. 

*) *r 

M » jet mach numbers * V./ c, 

1 1 


is the angle beyond which, theoretically , no sound is transmitted through the jet (4) 
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In the model development of the previous section, * nct is the angle at which the 
parameter T ^ goes to zero in equation 14. 

t 2 “ “ Hcos V 2 ' cos \] ** " 0 at *n " W 

where; 

ill ■ ~ - ft 

v n 2 

for the operating parameters in figure 3, the cut-off angle is: 


In theory, all the incident sound in the range >|i > ip ^ is transmitted through 
the jet, while no transmission occurs for sound incident at < i{> nct 

Refering to figure 3, the normalized sound pressure decreases more rapidly 
in the range of angles 40° < ij> n < 70° than it does on the range 70° < iji J1 < 90° . 
This indicates that transmission of noise through the jet is relatively efficient 
at <P n > 70°. while sound is transmitted through the jet at 4* n < 70°, the 
contribution of transmitted sound decreases rapidly. 

The normalized sound pressure contours reach a minimum at <f/ n ^ 40° . For 
angles less than 40°, the curves shift toward higher frequency „ In this range of 
angles, diffraction around the jet becomes dominant. The increase in normalized 
sound pressure with decrease in angle is consistent with barrier theory, where 
decreasing the angle of incidence has the same effect as increasing the distance 
between the source r-rt the barrier. 

The angle at which diffraction begins to dominate is evaluated from: 

i c „/ c. 

*ncd-“ s " ( T~ 3 


13 


This angle is the value at which the parameter in equation 14 goes to zero. 


M cos ■ 0 at tf) n ■ 


n v ncd 


where; 


if) * » 9 

r n 2 


for the operating parameters shown: 


♦ncd " 43,4 ° 

Since the transmitted noise decreases constantly for * n < Vt’ md diffracted 
noise adds to the sound pressure at angles * n < * ncd i the angle at which the 
noimalized sound pressure is minimum is equal to W 

Curves of normalized soimd pressure for the receiver in the transmission domi- 
nated, diffraction dominated find transition zones are shown in figure 4. The fre- 
quency scale is expanded to show the differences at higher frequency, 

For the receiver in the transmission dominated zone, the normalized sound 
pressure reduces rapidly, with frequency, to a minimum value. At higher frequencies, 
the sound pressure oscillates about a terminal value which is greater than the 
minimum. 

Mien the receiver is in the transition zone: 

^nct > ^n ^ncd 

the normalized sound pressure becomes smoother. The initial decrease in sound 
pressure is not as rapid. The first minimum occurs at lower frequency and is not 
as sharp. The sound pressure is oscillatory at higher frequencies, but the period 
of oscillation increases. 

For the receiver in the diffraction dominated zone, the normalized sound 


pressure reduction becomes more broad band. No oscillation at high frequency is 
observed. 

These general comments serve to identify the zones in which the mechanisms 
of transmission and diffraction control. The effects of jet operating parameters 
on the sound pressure in the transmission and diffraction dominated zones are dis- 
cussed below. 


2, Effect of Jet Operation Parameters 
A. Transmission Dominant Zone 

Figure 5 shows the noimalized sound pressure distribution plotted against 
non-dimensional frequency parameter, k Q a, for a receiver in the transmission 
dominant zone: 

900 > *n 5 *nct 
a) Effect of the Mach Number 

The curves on figure 5 represent a range of Mach Numbers, from 0.76 to 1.12. 
Generally, increasing the Mach Number shifts the curve toward lower frequency. 

Ai the Mach Number increases, the transmission zone cut+off angle , tj» t , is 
increased. Thus, a receiver at a fixed location is effectively moved closer to 
the cut-off angle; and, in so doing, sees less of the transmitted noise. The curves 
in figure 5 are found to collapse to one curve when plotted against: 

V C»n • 17 


\|i * observer angle (radians) 


H nct 


_ -1 , c o/ C;! 

cos f 

1+M 


h 


where : 


Figure 6 shows the resulting curve for the spacing ratio S/D ■ 2,667 from 
figure 5j as well as for spacing ratios, S/D - 1,000 and 5.000. 
b) Effect of Spacing Ratio 

As is shown in figure 6, the effect of increasing the spacing ratio is to 
shift the curves toward higher frequencies, These curves are found to collapse 
to a single curve when plotted against the normalized frequency parameter 

M <VW 17 £>" 17 

The resulting curve is shown in figure 7. This curve is similar to one 
obtained in a previous report (9) for the formulation of the two dimensional model. 
The onset of shielding occurs at the value of the normalized parameter approximately 
0.175 as was found previously. The minimum sound pressure occurs at a value of 
the parameter equal to 3.0, and was 2.72 in the previous analysis. The normalizing 
exponent, -.17, compares to the exponent, -.20, evaluated in the previous report. 

B. Diffraction Dominant Zone 

Figure 8 shows the normalized sound pressure as a function of frequency for 
a receiver in the diffraction dominant zone, 

^n < ^ ncd 

a) Effect of Mach Number. 

Figure 8 shows ‘chat the effect of increasing the Mach Number is to shift the 
curve toward higher frequencies. The curves approach a high frequency assymptotic 
curve which is independent of Mach Number. The curves in figure 8 are found to 
collapse to one curve when plotted against: 

-1.96 


k o a ^ncd - ♦n 3 


where : 


4> n ■ observer angle (radians) 

+ned ‘ diffraction zone “‘■off m S le 
-1 C n / C ' 

■ cos" ( — rr^-.) (radians) 

Mr 

C 

This relationship is valid only for M >• . 

c 

For M < — , \(i ^ does not exist; and there is no diffraction dominant zone. 

c * 

The curves in figure 8, for spacing ratio S/D * 2.667 are plotted against the 
normalizing parameter in figure 9. Also shown are similar curves for S/D ■ 1.000, 
and 5.000 . 

b) Effect of Spacing Ratio. 

From figure 9, the effect of increasing the spacing ratio is to shift the 
curves toward higher frequencies. As the distance between the source and the 
shielding jet increases, the minimum normalized pressure ratio which is approached 
at high frequency also increases. 

The minimum sound pressure ratio is found to be related to the spacing ratio 

% .59 log 1Q (S/D) + .15 
min 

c) Temperature 

i^ncd is related to jet temperature in that the speed of sound in the jet is a 
function of (T^K . Thus, as the jet temperature is decreased, the cut-off angle 
is also decreased. The expected result of decreasing the jet temperature is to 


uy; 


1L- 

Pin 
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shift the curve toward lower frequency. This result is shown in figure 10. The 
plots in figure 10 also show that the curves coalesce to a high frequency 
assymptate which is temperature independent. The assymptate is, however, dependent 
on spacing ratio. 

Figure 11 shows the two curves normalized against the parameter: 

V ( *ncd ' *n>' - ,96 

where: 

Vd - f CTJ.-S. 

The dashed line on the figure is the best fit curve from figure 9 for the 
corresponding spacing ratio. 

3. Comparison to Experimental Results, 

The analytical model is compared to experimental results by Kantola (10) 
and plotted in figure 12, The condition evaluated is : 

a. T- - 1238° R 

b. V-. * 1519 feet/sec 

J 

c. S/D * 2.667 

d. if> n * 30°, 60 c 

The analytical r»odel shows some agreement in form with the experimental data. 
The onset of shielding occurs at lower frequency for i|/ n * 30° than it does for 
\j» » 60°. As the jet axis is approached, the agreement between the analytical and 
experimental improves, with the significant difference being at high frequencies, 
The analytical model shows that normalized sound pressure increases for k Q a > 10. 
This indicates an upper frequency limit on shielding, The experimental results 
indicate no such upper limit in the range of frequencies investigated, 
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The major discrepancy between the analytical and experimental results occur 
at greater angles of incidence. This is felt to be due to the point source 
approximation of the source jet. The jet noise source is more realistically 
considered to be a distribution of point sources along the jet axis. The down- 
stream sources are contributing to the sound pressure measured in the shielded 
zone. From the previous discussion on transmission through the shielding jet, the 
sound incident on the jet at t|> < 70° is transmitted more readily than sound incident 

at \j» n < 70°. In general then, sound is transmitted upstream more efficiently than 
downstream because the angle of incidence is greater. Thus, the receiver at 

e ■ 60° "sees" more of the downstream noise sources than does the receiver 
T n 

located closer to the jet axis. The point source approximation for the source jet 
breaks down as the angle of incidence increases. 


V CONCLUSIONS 


In the shadow of the shielding jet in the downstream region, the mechanisms 
which control the noise are transmission and diffraction. 

Between the angles 90 > ^ > 'i' nct transmission through the jet dominates. 

In this zone, the normalized, sound pressure is close to 1. over a relatively large 
range of frequencies. As frequency is increased, the sound pressure decreases 
rapidly to a minimum value: at still higher frequencies, the sound pressure is 
oscillatory, indicating phase interference due to passage through the jet. The 
effect of increasing Mach or decreasing the spacing ratio number is to shift the 
curves toward lower frequency. The normalized sound pressure is represented by 
one curve when plotted against the normalized parameter 


(k 0 a) 



^nct-* x jj 


-.17 


The resulting curve compares favorably with a similar curve obtained from the 
two dimensional formulation. The value of the exponent, -.17, indicates that the 
sound pressure is relatively insensitive to changes in the jet operating parameters. 

Between the angles <ji t > ^ > ^ ^ , the influence of transmission decreases 
and the influence of diffraction increases. The normalized sound pressure curves 
shift toward lower frequency. The oscillatory nature of the sound pressure at high 
frequency is reduced. 

At angles, ^ < tji ^ , diffraction is the dominant mechanism. Since trans- 
mission constantly decreases for angles, ^ < ^ nc t 311(1 diffraction around the jet 

increases for >p n < , the minimum sound pressure occurs at tfs « ^ nc< j • The 

diffraction zone cut-off angle 


*ncd 


cos 
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is strongly dependent on the ratio of sound speed, , and the Mach Ntrher, M, 

For values, 

fo /c j > 1 
M 

^ncd cioes not ex i st » and the shielded region does not contain a diffraction dominant 
zone, 

In the diffraction dominated zone, the attenuation of the source noise is 
relatively broad band. At low frequency, the normalized sound pressure curves are 
shifted toward higher frequencies as the jet axis is approached. This result is 
consistent with barrier theory. At high frequencies, the normalized sound pressure 
approaches an assymptotic value, and none of the oscillatory character is observed. 
This indicates that at high frequencies, the sound incident on the shielding jet 
is scattered away from the shadow zone. 

In the diffraction dominant zone, low frequency reduction of the sound pressure 
is degraded as Mach Number is increased. At high frequency, the normalized sound 
pressure tends toward an assymptotic value which is independent of Mach Number. 
Increasing Mach Number increases the cut-off angle, 1J1 ^ . The shape of the 
normalized sound pressure curve is found to be related to the included angle be- 
tween the observer location and the cut-off, . When the sound pressure is 

plotted against the normalized frequency parameter: 

, „ , -,- 1.96 

V ( *ncd - *n ) 


the curves become Mach Number independent. 

Increasing jet temperature also increases the cut-off angle, since: 


c 


o 



c. 


3 



h 


Thus, when the normalized sound pressure is plotted against the parameter 

V «Wd • *„ >' 1,96 

the curves become independent of temperature, 

The value of the exponent, -1.96, indicates that the sound pressure in the 
diffraction dominant zone is more sensitive to changes in jet operating parameters, 
than is the sound pressure in the transmission dominant zone. 

The curves which result from normalization in terms of the cut-off angle are 
strongly dependent on the spacing ratio, S/D. The curves are shifted toward higher 
frequency as the spacing increases. This result is expected from barrier theory. 
The curves also appear to approach high frequency assymptotic values which increase 
as the spacing increases. The minimum sound pressure is found to be related to 
the spacing ratio by: 

% .59 log 10 ( S } + .15 

P in min 

for the range of spacing ratios 

1.000 < | < S.000 

The analytical model compares to experimental results favorably for angles of 
incidence close to the jet axis. The discrepencies that occur at greater angles of 
incidence are felt to arise from the difference between the single point source 
approximation in the model and the distributed source configuration of the real jet. 
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Figure 10. Variation of Normalized Sound Pressure with Non-dimensional Frequ 
Receiver in Diffraction Dominant Zone. Spacing Ratio, J| = 2.667. 






